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ABSTRACT 



Considered as the core of the Sco OBI association, the 
young open cluster NGC 6231 harbours a rich O-type 
star population. In 2001, the XMM-Newton satellite tar- 
geted the cluster for a nominal duration of about 180 ks. 
Thanks to the detector sensitivity, the EPIC cameras pro- 
vided an unprecedented X-ray view of NGC 6231, re- 
vealing about 600 point-like sources. In this contribu- 
tion, we review the main results that have been obtained 
thanks to this unprecedented data set. Concerning the O- 
type stars, we present the latest developments related to 
the so-called canonical Lx — iboi relation. The disper- 
sion around this relation might actually be much smaller 
than previously thought. In our data set, the sole mecha- 
nism that yields a significant deviation from this scheme 
is wind interaction. It is also the sole mechanism that 
induces a significant variation of the early-type star X- 
ray flux. In a second part of this contribution, we probe 
the properties of the optically faint X-ray sources. Most 
of them are believed to be low mass pre-main sequence 
stars. Their analysis provides direct insight into the star 
formation history of the cluster. 

Key words: Stars: fundamental parameters - Stars: 
early-type - X-rays: individuals: NGC 6231 - X-rays: 
stars - Open clusters and associations: individuals: NGC 
6231. 



1. INTRODUCTION 



Hot stars are known to be strong but soft X-ray emit- 
ters since the launch of the EINSTEIN satellite in Decem- 
ber 1978 (Harnden et al. 1979: Seward et al. 1979j. Al- 
though historically several hypotheses have been put for- 
ward to explain the origin of this X-ray emission, it is 
now commonly accepted that the latter is produced by 
shocks occurring within the denser layers of the wind. 
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These shocks, believed to result from the growth of insta- 
bilities of the line driving mechanism, heat the wind ma- 
terial to temperatures of the order of ten million Kelvin, 
thus producing a substantial amount of soft X-ray emis- 
sion. From t he observational point of view, it was soon 
realized (e.g. 'Harnden et aDIT979t iLong & White" 1980t 
Pallavicini et al. 1981) that the X-ray luminosities of 
hot stars were related to their bolometric luminosities 
through the so-called canonical Lx — iboi relation. In 
its generic form, this relation can be written as follows: 

Lx « IQ-'iboi. (1) 

The X-ray luminosities of O-type stars seem however 
to present a large scatter around this relation (e.g. 
iBerghofer et alJl997h . typically of about a factor 2 to 3. 

Beside the intrinsic emission of single stars, massive bi- 
naries are often m ore X-ray luminous compared to equiv- 
alent single stars (Pollock 1987: Chlebow ski & Garmaiivl 
1199 ll). This additional X-ray emission is usually inter- 
preted as the signature of a wind interaction resulting ei- 
ther from the collision of the winds of the two stars of 
the system or, if the wind from one component is much 
stronger than the other, from the interaction of the over- 
whelming wind with the secondary photosphere. This X- 
ray emission might further undergo phase-locked mod- 
ulations, resulting from the variation of the absorption 
along the line of sight. In eccentric binaries, such phase- 
locked modulations might also result from the change of 
the shock strength due to the varying distance between 
the two stars. 

Considered as the core of the Sco OBI association, the 
young open cluster NGC 623 1 hosts a rich early-type star 
population of different sub-spectral types and luminosity 
classes. Located at a distance of only about 1.6 kpc, it 
constitutes an ideal target to study the X-ray properties 
of an homogeneous (in terms of age, reddening, chemi- 
cal composition, ...) sample of OB-type stars. NGC 6231 
was thus chosen as the target of a 180 ks monitoring cam- 
paign performed in the framework of the Liege project for 
the Optical Monitor consortium guaranteed time. Car- 
ried out within a 5-day period in September 2001, the 
campaign actually consisted of six separate pointings of 
30 ks each. Nonetheless it forms one of the deepest X-ray 
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Figure 1. Combined EPIC-MOS three-color X-ray im- 
age of the young open cluster NGC 6231, showing the 
central part of the fov. It approximately corresponds to 
MOS CCD#1; the presented field is thus about 5' x 5'. 
The different false colors correspond to different energy 
ranges: red: 0.5-l.OkeV; green: 1.0-2.5 keV; blue: 2.5- 
lO.OkeV. The color image is also available from the 
XMM-Newton Image Gallery: http://xmm. vilsp a. esa. es/, 
external/xmm^cience/gallery/public/index.php 

views ever acquired of a young open cluster, but the par- 
ticular schedule of the campaign further allows to study 
the variability of the X-ray properties of the objects on 
different time scales. The analysis of t his remarkable data 
set is pr esented in a series of papers jSana et alJl20()5bl 
I2006alb^ to which we refer for further details. 

Beside the HMM-Newton observations, we have also 
monitored all the O-type stars and some of the bright- 
est B-type stars in the EPIC field of view (fov) by means 
of high resolution optical spectroscopy (mainly acquired 
using the ESO spectrograph FEROS). This data set has 
allowed us to provide further constraints on the ph ysi- 
cal pr operties of the early-type stars in NGC 623 1 ( Sanal 
120051) . In particular, we revised their multiplicity and re- 
derived their spectral classification. Thanks to the qual- 
ity of the FEROS data, we have detected the secondary 
signature for eight out of the nine O-type binaries in the 
EPIC fov. Together wit h the e xisting photometry of the 
cluster (e.g. Sung et al. Il998h . the tight constraints on 
the properties of the early-type stars provide a firm ba- 
sis for the X-ray analysis. The complementarity between 
the X-ray and optical data further constitutes one of the 
strengths of the present work. 



2. OBSERVATIONS AND DATA REDUCTION 

liMM-Newton performed the six successive exposures of 
our campaign during satellite revolutions 319 to 321. The 
fov was centered on the colliding wind binary HD 152248 



(•22000 = 16*^54'" 10!06, ^2000 = -41°49'30 '.'1), in the 
core of the cluster. Position angles (PAs) were very simi- 
lar through the six exposures. All three EPIC instruments 
were operated in the Full Frame mode together with the 
Thick Filter to reject UV/optical light. Due to the bright- 
ness of the cluster objects in the fov, the Optical Monitor 
was switched off throughout the campaign. 

The EPIC Observation Data Files (ODFs) were pro- 
cessed using the XMM-Science Analysis System (SAS) 
V 5.4.1 implemented on our computers in Liege. We ap- 
plied the emproc and epproc pipeline chains respectively 
to the MOS and pn raw data to generate proper event list 
files. No indication of pile-up was found in the data. We 
then only considered events with patterns 0-12 (resp. 0- 
4) for MOS (resp. pn) instruments and we applied the 
filtering criterion XMMEA_EM (resp. FLAG = 0) as rec- 
ommended by the Science Operation Centre (SOC) tech- 
nical note XMM-PS-TN-43 v3.0. For each pointing, we 
rejected periods affected by soft proton flares. For this 
purpose, we built light curves at energies above lOkeV' 
and discarded high background observing periods on the 
basis of an empirically derived threshold (adopted as 0.2 
and 1.0 cnts^^ for the MOS and pn instruments respec- 
tively). The so-defined GTIs (Good Time Intervals) were 
used to produce adequate X-ray event lists for each point- 
ing from which we extracted images using x- and y-image 
bin sizes of 50 virtual pixels ^. 

We finally combined the event lists obtained for all six 
pointings to improve the statistics of faint sources. For 
this purpose, we used the SAS task merge. For each 
EPIC instrument, we included the event lists resulting 
from different pointings one by one. We also built merged 
event lists that combine the twelve MOS or the eighteen 
EPIC event lists. The Attitude Files generated by the 
pipeline were merged using the same approach and we 
adopted, for handling the merged event lists, the Calibra- 
tion Index File (CIF) and the ODF corresponding to the 
first pointing. 

The total effective exposure times towards the cluster 
are, respectively for the MOSl, MOS2 and pn instru- 
ments, of 176.5, 175.0 and 147.5 ks. Together with the 
high sensitivity of the XMM-A^ewfo/i observatory, the 
combination of the six pointings and of the three in- 
struments provides one of the deepest X-ray views of a 
young open cluster Fig. ^ shows a three-colour image 
of NGC 623 1 and reveals a densely populated field with 
hundreds of point-like X-ray sources. We estimate our 
detection flux limit to lie between about 3 x 10"^^ and 
L5 X lO^^^ergcm^^ s^^ depending on the position on 
the detectors and on the source spectrum. 

For each source, we finally adopted a circular extraction 
region with a radius corresponding to half the distance 

'Expressed in Pulse Invaiiant (PI) channel numbers and consider- 
ing that 1 PI channel approximately con'esponds to 1 eV, the adopted 
criterion is actually PI > 10000. 

^Though the physical pixels of the EPIC MOS and pn detectors 
have an extent on the sky of respectively 1'.' 1 and 4'.' 1, the virtual pixels 
of the three instruments correspond to an extent O'.'OS. The obtained 
images have thus a pixel size of 2'.' 5. 
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Figure 2. ISM-absorption corrected X-ray luminosities in 
the 0.5-lO.OkeV band plottedvs. bolometric luminosities. 
The dijferent symbols indicate the dijferent properties of 
the sources. Spectral type: O (filled symbols), B (open 
symbols). Luminosity class: supergiant (black), giant 
(red), main sequence (blue). Multiplicity: binary (trian- 
gles), presumably single RV-variable star (squares), pre- 
sumably single RV-constant star (circles). Best-fit linear 
relations in the log — log plane for O (Eq. |2j and B (Eq. 
stars are indicated by the dashed and the dotted lines 
respectively. The two vertical bars in the upper part of 
the graph give the expected 1-a deviation for B (a^) ond 
O ((To) type stars. 



to the nearest neighbouring X-ray source. Due to the 
crowded nature of the cluster core in the X-rays and to the 
limited spatial resolution of the EPIC detectors, the back- 
ground could not be evaluated in the immediate vicinity 
of the stars, but had to be taken from the very few source 
free regions. The details are given in S ana et al. (2006ai) . 
We then used the appropriate redistribution matrix files 
(rmf) provided by the Science Operations Centre (SOC), 
according to the position of the considered source on the 
detectors. We built the corresponding ancillary response 
files (arf) using the ARFGEN command of the SAS soft- 
ware. The spectra were finally binned to have at least 10 
counts per bin. 



3. THE OB STAR POPULATION 



Prior to our scientific analysis, we first went through a 
deep overview of the existing literature and databases on 
NGC 623 1 and we made a census of the early-type star 
population in the observed fov. Our census mainl y relies 
on selected published works (see references in Sana et al. 
l2005 b). on the Catalog of Galactic OB Stars ( Reed.2003.) 
and on the WEBDA^ and SIMBAD"* databases. This re- 



sulted in 108 objects among which 92 B-stars, 15 O-stars 
and one WR system (WR 79). 

All the O-type stars were detected by XMM-Newton as 
soft but strong sources, which allowed us to study the 
complete population rather than a subsample of it. On 
the other hand, only about 20% of the B-type stars could 
be associated with an X-ray source. Using the XSPEC 
software, we adjusted the X-ray spectra of the differ- 
ent sources associated with early-type objects. In doing 
so, we used up to three-component thin thermal plasma 
models (mekal models) allowing for a distinct local ab- 
sorption column (wabs model) for each component and 
accounting for additional absorption by the interstellar 
medium (ISM). The latter column was held fixed to a 
value computed from the reddening of the different ob- 
jects, using the typical colo urs o f Schmidt-Kaler ( 1982) 
and the gas to dust ratio of Bohlin et alj (il978 ). Using 
the best-fit models, we finally computed the X-ray lumi- 
nosities, corrected for the ISM absorption column only. 
We also recomputed the bolometric luminosities, adopt- 
ing our revised sp ectral classification an d the bolometric 
correction scale of Schmidt-Kaled i 1 982^ . Fig.|2]presents 
the location of the different objects in the Lx — iboi di- 
agram. Note that Fig.|2]is restricted to objects fitted with 
at least 2-tempe rature mekal com ponents. Reasons for 
this are given in lSana et'al] (|20063). 



3.1. The O-type stars 

Focusing in a first step on the O-type stars, we note a 
clear linear relation in the log — log plane. However sev- 
eral points deviate significantly from th is canon ical rela- 
tion. HD 152248 (07.5111(f) + 07III(D. ISanaet al. 2001.) 
presents clear evidence for a wind-wind interaction that 
we traced both from the optical and X-ray domains. 2-D 
hydrodynamical simulations of the collision further rea- 
sonably reproduce the observed modulations of the X- 
ray flux of the system, l ending further support to this in- 
terpretation (Sana et al j|2004i) . CPD-4r7742 (09V + 
B1.5V, Sana et al. 2003) X-ray light curve also displays 
clear though apparently complex modulations of its X-ray 
flux, that we have however related to a wind interaction 
of a peculiar kind. In this Oh-B system, the dominant pri- 
mary wind most probably crashes into the secondary sur- 
face. A simple phenomenological model of such a wind- 
photosphere interaction indeed reproduces pretty well the 
main features of the X-ray light curve (Sana et al. 2005a). 
It is thus clear that the two latter systems are not repre- 
sentative of the intrinsic X-ray emission of single O-type 
stars. Excluding these two points, a least-square linear fit 
yields; 



logix - logiboi = -6.912(±0.153) 



(2) 



^httpV/obswww.unige.ch/webda/ 
''http://simbad.u-strasbg.fr/Simbad/ 



where the X-ray luminosity is given in the 0.5-10.0 keV 
band. It is interesting to note that we also adjusted a 
power-law relation (thus in the form logix = F x 
log Lhni + K). How ever, as indicated by a i^^-test (see 
e.g. lBevingto'nlll969h . this additional parameter does not 
improve significantly the quality of the fit compared to 
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Figure 3. Hertzsprung-Russell diagram o f the EPIC 
sources with optical counterparts in the 'Suns e t al\ 
{J998) catalogue. Evolutionary tracks from Siess e t al\ 
hoOO ) for masses of 0.2, 0.3, 0.4, 0.5, 0.7, 1.0, 1.5, 2.0, 
2.5, 3.0, 4.0, 5.0, 6.0 and J.OMq are overplotted. Filled 
dots, filled diamonds and open triangles indicate respec- 
tively Ha emitting stars. Ha candidates and stars with 
no evidence for emission. The thick solid line shows the 
ZAMS, while the dashed lines correspond to isochrones 
forages of 0.5, 1.5, 4.0, 10.0 and 20.0 My r 



the scaling law of Eq. |2l In the following, we thus de- 
cide to adopt the simplest form of the canonical rela- 
tion, as quoted in Eq.|2] One of the most outstanding re- 
sults lies however in the limited dispersion (of only about 
40%) of the X-ray luminosities around the best fit rela- 
tion. Compared to the scatter (a factor of 2 to 3) observed 
by^erghofer et al. ( 1997), this suggests that the intrinsic 
X-ray emission of 0-type stars might be much more con- 
strained by the physical properties of the star than pre- 
viously assumed. It is also remarkable that, within our 
sample, only the binaries present a significant modula- 
tion of their X-ray flux, suggesting thus a wind interac- 
tion origin. Excluding the case of HD 326329, proba- 
bly contaminated by a neighbouring flaring source, extra 
emission produced in a wind interaction region is also the 
only mechanism that yields a significant deviation from 
the canonical Lx — iboi relation. All in all, our analysis 
suggests thus that the intrinsic X-ray emission from sin- 
gle O-type stars might be rather stable, both in terms of 
time-variability and of scattering compared to the mean 
Lx — iboi relation. 
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3.2. The B-type stars 



Figure 4. Distribution of the ages of X-ray selected PMS 
objects as interpolated from the isochrones derived from 
the Siess et al. (2000) evolutionary tracks. PMS candi- 
dates with A(i? — Ha) > 0.12 are indicated by the solid 
line, while those with A{R — Ha) < 0.12 are indicated 
by the dashed line. The total numbers of objects with 
A{R-Ha) > 0.12 and A{R- Ha) < 0.12 are respec- 
tively 93 and 303. 



dence level of 0.99 in favour of the presence of a correla- 
tion between logLx and logLboi- A linear least-square 
fit yields: 

logLx = 0.22(±0.06) logLboi + 22.8(±2.4) (3) 

However, we emphasize that the undetected B-type stars 
were not taken into account in the derivation of the 
present relation. The fact that about 80% of them lie 
below our detection threshold suggests that either X-ray 
emission from B-type stars is not an intrinsic property of 
such stars or, at least, that it is not fully governed by their 
bolometric luminosities. In the latter case, we could only 
be detecting the upper envelope of a largely scattered dis- 
tribution in such a way that, by coincidence, the detected 
sources suggest the presence of a linear relation. Beside 
these considerations, one may further note that about one 
third of the detected B sources present flaring-like activi- 
ties during the time-span of our XMM-Newton campaign. 
Their X-ray spectral properties are further reminiscent of 
those of PMS stars. It is therefore difficult to definitely 
conclude whether the detected X-ray emission is directly 
associated with the B-type stars or if it is rather produced 
by PMS stars either as the secondary component in a bi- 
nary system or being located by chance along the same 
line of sight. 



Turning to B-type stars, their distribution in the log Lx ~ 
log Lboi plane also suggests the presence of a linear rela- 
tion, although as quoted above, only about 20% of these 
objects are seen in the X-rays. Indeed the linear corre- 
lation coefficient is r ^ 0.75, corresponding to a confi- 



4. THE OPTICALLY FAINT X-RAY SOURCES 



Beside the early-type stars, the EPIC images reveal sev- 
eral hundreds of additional point-like sources that clus- 
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Figure 5. Locations of the NGC 6231 0-type stars in 
the H-R diagram. Symbol shapes and colours provide 
different indications on the object nature: black, super- 
giant; red, giant; blue, main sequence; triangle, be- 
longing to a binary system (filled triangles indicate the 
heaviest star of the system while open triangles stand for 
the less massive companion); square, presumably single, 
RV-variable star; circle, presumab ly single, constant-RV 
star The evolutionary tracks from \Schaller et al.\ il99^) 
have been plotted (plain lines) together with isochrones 
(dotted lines) computed for ages ranging from 2 to 10 
Myr with a step of 2 Myr. The figure appears in colour in 
the electronic version of this work. 



ter towards the core of NGC 6231, suggesting a phys- 
ical link with the cluster. T hese are typically fainter but 
harder than the O-type stars JSana et al.l2005bl) . As a first 
step, we compared the X-ray source list with different 
existing optical/lR catalogues, adopting a limited cross- 
correlation radius of 3". As a result, almost 80% of the X- 
ray sources in the XMM-Newton fov could be associated 
with at least one optical/lR counterpart. In particular, we 
made use of an extended version of the UBV(RI)c Ha 
photometric catalogue of Suns et al. ( 1998). Roughly, 
the latter work covered a square field of 20' x 20' that is 
inscribed within the 15' radius fov of the EPIC cameras. 
NGC 6231 being a quite young open cluster (age^3- 
5 Mvr. lBaume et al.l 19991 ). we expect that the lower mass 
stars have not reached the ZAMS yet. Being probably 
still embedded in their natal clouds, they might suffer a 
heavy circumstellar absorption. X-ray emission is how- 
ever one of the best s election criteria for these objects 
dSung & Besselll 12004*) and the present XMM-Newton 
image of NGC 623 1 might form one of the best censuses 
of the cluster PMS population so far. Ha emission is a 
further well known property of the classical T Tauri stars. 
Among the 53 6 X-ray sources lo cated within the field 
covered by the ISung et al] ( Il998h catalogue, we indeed 
identified 93 candidates displaying a significant excess in 
the R—Ha col our index (adopted as A(R—Ha) > 0.12, 
see ISana et alJi2006h. for details). 



Fig. |3l presents the location of these optically faint X- 
ray sources in the Hertzsprung-Russell diagram. Most of 
them are indeed located above the ZAMS as expected for 
PMS objects. Comparing with PMS evolutionary tracks 
of iSiess et alJ (j2000j), it appears that most of them are 
low-mass stars (M < 2 Mq) that started their forma- 
tion about 2 to 4 Myr ago. The distribution of the ages 
of these X-ray selected PMS stars is presented in Fig.|3 
and suggests that star formation in NGC 623 1 might have 
started at least 10 Myr ago at a relatively slow rate. The 
latter then slowly increased to culminate in a starburst- 
like event about 2 to 4 Myr ago, a period that also corre- 
sponds to the formation of the cluster massive stars (see 
Fig. |5). It is further interesting to note that, neither in 
Figs. |3]nor0] a clear difference appears between the Ha 
emitting X-ray sources (probably classical T Tauri stars) 
and those that do not present evidence of such emission 
(probably weak-line T Tauri stars). No difference could 
also be found in their respective spatial distribution in the 
studied fov. 



5. CONCLUSIONS 



We briefly presented some of the main achievements 
of an XMM-Newton monitoring campaign of the young 
open cluster NGC 623 1 in the Sco OB 1 association. Of 
a nominal duration of 180 ks, it was actually split into 6 
separate pointings spread over 5 days. Clearly it consti- 
tutes one of the deepest X-ray observations of a young 
open cluster and its particular scheduling allowed us to 
probe the variability of the detected sources on different 
time scales. The EPIC cameras revealed a crowded fov 
with more than 600 X-ray sources. The large majority 
of these could be identified with an optical/IR counter- 
part and their location in the H-R diagram suggests that 
most of them are actually low-mass (M < 2 M0) PMS 
stars. Their study seems to indicate that the star formation 
in NGC 6231 was probably not an instantaneous event 
but might have started at least 10 Myr ago at a relatively 
slow rate. The bulk of the cluster stellar population has 
however started its formation during a starburst-like event 
about 2 to 4 Myr ago, an epoch during which the most 
massive stars of the cluster were also formed. 

XMM-A^ewfo« also detected the complete O-type star 
population in the EPIC fov. Being strong but soft emit- 
ters, these objects clearly dominate the EPIC images. 
Restraining our analysis to the 0.5-10.0 keV band, we de- 
rived a new value for the Lx — ^boi canonical relation, 
expressed in the form of a scaling law. One of the most 
outstanding results is the limited dispersion of the O-type 
star X-ray luminosities around this new canonical rela- 
tion compared to the one observed by previous studies. It 
is clear that the fact that we have been able, among oth- 
ers, to identify and reject the colliding wind binary sys- 
tems from the fit has helped to reduce the scattering. The 
particular homogeneity of our sample, notably in terms 
of age, metallicity and reddening, might have played a 
crucial role in this regard. 
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Finally, it is also worth to note that, within our sample, 
variability of the X-ray flux is only observed for probable 
colliding wind binaries. It is further the only mechanism 
that seems to produce strong deviations from the canoni- 
cal relation. All in all, our analysis suggests thus that the 
intrinsic X-ray emission from single O-type stars might 
be much more stable than previously thought, both in 
terms of time- variability and of deviation from the mean 
Lx — iboi relation. 

The comparison of the present results with others, derived 
in an homogeneous way from the study of other young 
open clusters, will therefore be of particular interest to 
probe the universality of the canonical relation. It could 
further help to investigate the influence of other important 
parameters, such as age and metallicity, which might cru- 
cially affect the intrinsic X-ray emission of O-type stars. 

Finally we note that, compared to XMM-Newton, future 
X-ray missions should definitely combine a high sensi- 
tivity with an increased spatial resolution. The latter is 
indeed crucial in very crowded fields such as the one stud- 
ied in this contribution. We also note that, while deahng 
with such an extended data set, 12 months is a very short 
time. One can thus wonder whether the proprietary pol- 
icy could not be extended to allow scientists to perform 
in-depth studies in good conditions. 
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